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ABSTRACT 


A method is presented for the analytic design of linear 
control systems. Design is accomplished by the solution of 
an algebraic set of simultaneous equations derived from the 
design specifications. Specifications considered include the 
system dam#iing ratio, undamped natural frequency, bandwidth, 


steady state error coefficients and root sensitivity. 





The method is shownste areal, accurate, specific and 
capable of working with any number of parameters. In 
addition, the analytic approach minimizes dependence on 


graphical techniques and trial and error solutions. 
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SYMBOLS AND ABBREVIATIONS 


Mow 


Superscript or 
Subscript "o" 


w (omega) 


Wn 


Characteristic equation 
Steady state error 

Forward transfer function 
Acceleration error constant 
Accelerometer gain 

Gain 

Position error constant 
Velocity error constant 
Tachometer gain 
Reasonance peak 


Specified or numerical value 


Angular frequency 

Undamped natural frequency 
Bandwidth ( -3 db) 

Reasonant frequency 

Polynomial in & 

General sensitivity parameter 
Complex variable of the form ao + jw 
Sensitivity 


Sensitivity of P with respect to Q 


¥F component of root sensitivity 


wy, component of root sensitivity 
Total root sensitivity 


Damping ratio 





1. Introduction: 

Linear control system design can be defined as "a set of 
procedures for determining the value of system parameters 
which will cause a physical system to operate in a prescribed 
manner." The set of procedures noted above is as broad as 
the definition. One may choose to perform the design effort in 
the time domain by applying modern state space concepts or by 
direct solution of the system differential equations. Alternatively, 
the complex frequency domain may be selected, employing any 
of the universally accepted procedures such as Nichols and 
Bode plots, root locus, Mitrovic's Method, etc. 

Notwithstanding the achievements which have resulted 
from direct application of the above procedures, it must be 
admitted that they have inherent limitations. With the 
exception of simple systems, time domain analysis virtually 
demands the assistance of a digital computer. Ina sense, 
this is a disadvantage because computers are neitler available 
to all people at all times nor is it economically practical to 
employ a computer for all computation. The complex frequency 
schemes are considered imperfect because they are time 
consuming, i.e., essentially graphical in nature and often 
based on trial and error, and normally limited to one or two 


variables, 


In view of the above discussion, it appears that there is 
need for still another design tool. However, in order fora 
new procedure to even merit consideration, it must be 
relatively free from the faults and limitations cited above. 
Specifically, the method should be analytically simple, i.e. 
extend the range of problems which can be solved without 
resorting to graphical techniques or the exclusive use of a 
computer. It should also be rapid, capable of handling several 
variables, specific and accurate. The remainder of this paper 
is devoted to the development and discussion of a design 


procedure which is believed to have these attributes. 





2. The Algebraic Method: Concept and Development 


Concept 

The proposed design procedure is based on the premise that 
system design specifications can be transformed into algebraic 
functions of the system parameters. If the specifications are 
independent, N specifications will produce N equations which 
simultaneously constrain the system. Control system design 
would thus revert to the relatively simple problem of solving a 
set of simultaneous algebraic equations. 

The problem of transforming specifications into equations 
is not as broad and difficult as one might initially assume. 
The most common and practical complex frequency specifications 
are: undamped natural frequency (w,), daming ratio (—), steady 
state error €..) , bandwidth (wy) , reasonance peak (Mow) and 
reasonant frequency (w,). In many cases, both ¥ and w, are 
specified. This is equivalent to specifying the location of a 
complex pole pair and generally implies that these poles should 
dominate the system response. If attention is restricted to this 
situation, then the specifications of primary interest are ; 
Wy, Wp and Eos because Mow and w, can be at least 
approximately expressed in terms of © and Wh: The remainder 


of this section is therefore devoted to the development of methods 


which can be used to obtain system constraint equations from 
specifications for F , wy, wy and Egg. 
Root Location (¢ and w,) 

The coefficient equations associated with the specification 
of a complex pole pair can be developed in two ways. The first 
method is based on the root~-coefficient relations of an 
algebraic equation; the second is based on Mitrovic's Method. 

Assume the systems' characteristic equation (CE), obtained 
by setting the denominator of the closed loop transfer function 


equal to zero, is the polynomial. 


v1 
= An Sk =O 
=O 


(2-1) 
where s is a complex variable and defined as 
S =— FW, + 4 Wnfi— EX (2-2) 


and the as represent the real coefficients which are functions 
of the system parameters. 
Since equation (2-1) is of order n, there must be n roots 
which satisfy the equation. Therefore, an equivalent expression 
s. 
for equation (21) is 
7 
(s + Ri) = © 
¢= / (2-9) 


where R,; represents the ith root of the polynomial. By expanding 


equation (2-3) and equating the coefficients of like powers of s in 


1c 





equations (2-1) and (2-3), n equations can be obtained which 
relate the roots and coefficient of the original polynomial. 
If & ‘ and wo define the location of the specified complex 
pole pair, it is permissible to set 
Ki = foWo +f wovl— For 
(2-4) 
R2 = Fowe — j We Vi— we 
(2-5) 
Substituting these values of Ry and Ry into the n root~coefficient 
relations obtained previously gives n equations in terms of the 
remaining n-2 roots and the coefficients of equation (2-1). By 
algebraic manipulation, it is then possible to successively 
eliminate each of the remaining n-2 roots and finally obtain 
two equations involving only the coefficients. Since the 
coefficients are explicit functions of the system parameters, 
any combination of the parameters which satisfy these equations 
will generate the desired pole pair. 
Example 2.1 


Let the system characteristic equation be 


SF tars +Q:S + Qe =O 


Find the relations which must exist between the coefficients in 
order to locate a pole pair with = 1/2 and wo = 10. 


Denote the three roots of the CE as Rj, Rg and R3. Then 
(s+®,)(s+R2)(St+Rs) = $s? +925*4+4%5 +0 
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and a2, =R,R2 + Ra (Ri + R2) 
Geo = KF, K2 Ke 
Gz = Ki rR2 +Kz2 


Setting Ri = & + J/0 J Yu. 
R2 = 5S- ji0 (Hd 

Then /20 (Gz-/0) =Ge (1) 
10(Qi-/0) =Qe (2) 


/O0Fs = Go 





Equations (1) and (2) are the desired coefficient relations; 
equation (3) provides information about the third pole. 

While the root-coefficient method works well with low 
ordered systems (and even provides information about the 
remaining n-2 roots), it requires an inordinate amount of 
computation when applied to systems of higher order because 
the coefficient relations must be extracted from a set of non- 
linear equations. This difficulty can be overcome (at the expense 
of losing information about the remaining n-2 roots) by using 
Mitrovic's Method. With this method, a kth ordered system can 
be handled as easily as one of third order. The method is fully 
developed in (3) and (5). Only a brief review is given here. 

The argument begins with equations (2~1) and (2-2) defined 
previously. Two independent equations in a,, Wy and fircan be 
obtained by substituting equation (2-2) into equation (2-1) and 


setting the real and imaginary portions of the resulting 


V2 





expressions equal to zero. Through additional manipulation and 


the introduction of a new variable, the equations take the form 


z Qn Wn bk (F)=0 (2-6) 
K 
”") oy =O 
= Qe a q :. ae -7) 
where p, (F) =— [2F 4, (F) +f, (F) 7 (2-8) 


with @ (F)=0, A(F)=-1 ond 4, (p= -4n CF). 
Other values of Dr, ( F ) for various values of ¥ are tabulated 
in Appendix I. 

Equations (2-6) and (2-7) are known as the generalized 
Mitrovic equations and can be used to obtain the coefficient 
relations required to locate a pole pair as specified by F A 
and w,. 

Example 2.2 
Work example 2.1 using Mitrovic's Method. 


For a third order system, the generalized Mitrovic equations 


are 
” 2 3 
Heh (F) +4,0, 08) + a4, (F)w, + 9,0, £(F) =O 


208, (F) +a, wd (F) +92 we p, (F) + a3 WF, $26 F)= O 


Substituting wp = Wo = 10 and using values for 4,(1/2) from 


Appendix I, the above equations simplify to 


= 


Oe = /o Go=/0002-/000 


The results are equivalent to those obtained previously. 
Although only a single pair of complex poles is being 
considered here, it should be noted that either method is 
capable of handling any number of real or complex poles 
(providing, of course, that the total number considered is less 
or equal than the order of the system). For example, assume 
two complex poles pairs and two real roots of a sixth order 
system are specified. The root-coefficient method will provide 
six constraining equations by direct substitution of the given 


values. If Mitrovic's Method is used, four constraint equations 





dependent on the complex poles can be obtained by repeated 
use of equations (2-6) and (2-7). The remaining two equations 
which are dependent on the real roots are produced by repeated 
use of 
“ K 
= Cl) an GR = O 
K=0 
(2-8) 

obtained by letting s = ~g-in equation (2-1). 

Steady State Error (Eee) 
There are three different steady state error constants which 


may be specified. They are 


position error constant = Kp = lim G(s) (2-9) 
SO 
velocity error constant = mo = lim sG(s) (2-10) 
Seo 
acceleration error constant =k, = lim s2G(s) (2-11) 
Sa-O 
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The above definitions are only valid for systems which have all 
their roots in the left~half of the s~plane (because the 
definitions are merely restatement of the final value theorem 
for complex variables) and unity feedback around the forward 
loop transfer function, G(s). To satisfy the last requirement, 
it may be necessary to manipulate the block diagram into an 
equivalent G(s) block with the required unity feedback. 

The constraint equation is obtained by equating the specified 
value of Eee with the algebraic expression resulting from 
application of equation (2-9), (2-10) or (2-11). The procedure 
is demonstrated in the following examples. 

Example 2.3 

For the system shown in figure 2-1, obtain an expression 
between K, Z and P for Ss = 6, 

For the given system 

Gee (S)-0=. ne K (S+ @) 
s2 (s+5)(5+6)(S+*P) 


Applying equation (2-11) 


= ee “- A hey 
#, —" Gea) 30P 


which, for Ea o, Teguceci to: 180P = KZ 
Example 2.4 
For the system shown in figure 2-2, derive an equation in 


terms of Ky. Ko and K, which will give Ky = 8, 
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Since the system does not have unity feedback, insert 
fictitious positive and negative feedback paths around the 
system. Geq(s) can then be found by including the positive 


feedback .branch. 


Ky 


Gee CS) = 
a St(1+K,Kk3) +38 (S4+K, Kz) —KK, 


Applying equation (2-10) 
= L — 
Ky £iyn s Gea(s) = O 


Therefore the specification cannot be satisfied. 
Bandwidth (wy) 

Bandwidth will be defined as the angular frequency, Wp: 
at which the magnitude of the closed loop transfer function has 
a value of 0.707, i.e., - 3 db bandwidth. The transformation 
of the specification into the desired constraint equation can be 
accomplished by following the usual procedures: 

A. Find the closed loop transfer function and set s = jw. 

B, Find the magnitude of the resulting expression and 
equate it to 0.707. 

C. Substitute w = wy and simplify. 

Example 2.5 

For the system shown in figure 2-3, derive a relation in 

terms of K and K, which assures the system has a bandwidth 


of 10 radians/second. 
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Following the prescribed procedure: 





~<S)_ = K 
Rs) S24 Ss +SKKr tK 
Cyw) F_ Kz 





Rgyw) | (K-w? 2 + w2(s+KKr)2 
Por w = wy = 10, [eee |° , Wi 
KR (jw) 


o Ke = s0t—200K +/0O0(S+KEr)* 


The constraint equation obtained in the preceeding example 
for a second order system is typical of the expressions which 
are developed from bandwidth specifications. Since the 
expressions increase in complexity for higher ordered systems 
or for systems having closed loop zeroes, most design problems 
involving bandwidth specifications will usually require the 
assistance of a computer or a graphical technique. (An excellent 
graphical procedure based on Mitrovic's Method is given in (4) .) 
If such aids are not available, a convenient, and relatively fast, 
alternative solution is to initially ignore the bandwidth 
specification and design the system to meet the remaining 
specifications. The results can then be inserted into the system 
bandwidth equation to determine if they also satisfy the bandwidth 
requirement. if the specification is not satisfied, the magnitude 
of the deficiency may suggest a more reasonable structure for the 


next attempt. 


l/ 
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BLOCK DIAGRAM FOR EXAMPLE 2.3 
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[ K3s2 _ 
FIGURE 2.2 


BLOCK DIAGRAM FOR EXAMPLE 2.4 
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PIGURE 2.3 


BLOCK DIAGRAM FOR EXAMPLE 2.5 
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3. Application of the Algebraic Method. 

As stated previously, the basis of the proposed design 
method is the transformation of N specifications into a set 
of N simultaneous constraint equations that are explicit functions 
of the system parameters. Since the specifications are 
independent, a specific design will generally require that the 
system have N parameters. The insertion of sufficient 
parameters, usually compensating Bienes. is an integral part 
of the design procedure. 

It should be noted that the mere selection of sufficient 
parameters does not guarantee a satisfactory solution. A 
converse of this statement is also true. An unsatisfactory 
solution does not imply that a solution does not exist; one may 
possibly be obtained by employing a different set of parameters, 
These statements suggest that there is some trial and error 
associated with the procedure. This is correct. However, the 
amount of trial and error required can be minimized through 
experience and application of control theory fundamentals. In 
this respect, root locus sketches are particularly valuable in 
obtaining clues to the nature of the required parameters. 

The algebraic method has the inherent capability to handle 
any number of specifications or parameters. Therefore a 
parametric solution can generally be obtained by simply inserting 


more parameters into the system than there are design 
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specifications. If M parameters are selected for use in the 
design of a system subject to N specifications, for M>N, the 
results will be parametric and expressable in terms of any M-N 
parameters. 

Example 3.1 

Given G(s) = K/s(s + 4), design* a system which will 
produce a pair of complex poles defined by F , = 1/2 and wy = 4 
wilt Kk. =f 

It will be stated without proof that K = 16 satisfies all 
specifications. 

This seemingly trivial example has major significance. 
Specifically, one parameter successfully satisfies three 
specifications which is a direct contradiction of the previous 
discussion. It would appear that adherence to the general rule 
of N parameters for N specifications could possibly result in 
the incorporation of redundant material into the system. This 
does not seem likely because of the following observation: 

In all cases studied, parameters which did not assist in 

the solution tended to be eliminated. 

Examples 3.2 and 3.3 will be used as illustrations of this 


observation. 





Example 3.2 
Insert cascade compensation, G(s) = (st+Z)/(s+P), into 


the system given in example 3.1 and redesign the system to 


*Note: All designs will employ unity feedback between output and 
input (in addition to any other specified feedback). 
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meet Lhe sene srecii le *vions. 


S3 + S27(4¢+P) + 4Ps 


Applying equation (2-10): 


a z (1) 
Ky = 
System CE: 
St+st(4+P) +5 4P+K) + Ke =O 
Mitrovic equations: {(4+P) = yYyPrK (2) 


KZ = /6(4+P)-64 (3) 

A simultaneous solution of equations (1), (2) and (3) gives: 

K= 16 Y =aa 
Since the pole and zero of the filter are coincident, the filter 
performs no function and can therefore be eliminated. 

Example 3.3 

For G(s) = 50/s(s + 5), it is possible to locate poles at 
-52j5 (PF, — 0. M7 and w, = 5 (2) 172) by using either 
tachometer feedback (K, = 0.1) or a cascade compensator (Z = 5, 
P= 10). In each case, K,,=5. Design a system using both 
tachometer feedback and a series compensator to satisfy the 
same specifications. 

For combined tachometer feedback and cascade compensation: 


SO (St+#) 


3) = aaa" 
Gee S3 +S2(P4S54S0Kr) +5(5Pr S0KT#) 


Applying equation (2-10): Ky = s0z/(5P +50 Kr2) (1) 
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System CE: : 


S33 4+ 52 (Pr57+50Kr)+ 5S (SP+50Kr2 +50) + S02 =O 


Mitrovic equations: 402 = 50 (P+ 5% S0Kr)-S00 (2) 


SP rS50Kk72 +505 /0(P+E5 +r SO0Ky;)-50 (3) 


Simplifying: z = 5 (4) 
P=10(1 - oK,) (5) 

The result obtained in this example is doubly significant. 
First of all, it indicates that the procedure is capable of 
detecting conditional as well as complete redundancy. In 
example 3.2, the filter was redundant for all values of the 
compensator pole and zero. In this example, however, 
equation (5) shows that the filter is redundant only when 
K,= 0.1. Similarly, the tachometer may be discarded when P = 10. 
The second point of interest is the fact that a unique solution 
may not exist. In this case, any positive values of P and kK, 
which satisfy equation (5) also satisfy the specifications. Which 
solution is best? At this point, it is a matter of conjecture. 
However, this example will be rediscussed in Section 4 after a 
means for determining the best design has been developed. 


Example 3.4 


K 


Given the forward loop transfer function, G(s) = Ss+ in) Gee J 


design series compensation to give roots with F , = 1/2 and 


Wo = 6. The minimum acceptable value of Keres 4.15. 


ae 











For cascade compensation: 


K (s+2#) 
—) 
Gee oe s*+7S5P 334 S2(75P +3.5) + 3eS PS 
Applying equation (2-10): Kyv=4./5 = se (1) 
3-5P 
System CE: 


ste 75 Ps3 + 5? (25 P4385) +5(3-5 PK) + KZ=0 


Mitrovic equations: KE = SY P-/49S (2) 


K = 4AILSP- 19S (3) 


Solving equations (1), (2) and (3) gives: K= 1475 


Z= 0.373 
Example °3.5 
The bandwidth for the system in example 3.4 is about 6.5 
radians/second. Redesign the system so that the bandwidth is 
reduced to 4 radians/second. 
The introduction of another specification requires the 
insertion of another parameter. Use tachometer feedback around 


the original plant and series compensator. Then 


= K(S+#2) 
Geq(s) a i oa ——— ss =. 
S++ 53 (7, 570) + 5% (3e517.5P KK) + S$(BSP+KKr=) 
Applying equation (2-10): Kv =“4.1S = 2] (1) 
3S P+KKrz 


System CE: 
sty+ 59 (Z7S54rP) + 8% (3-5$4+7.5P4KKr) +s (305P+KKr®rk) +K2E=0 


Mitrovic equations: ps 2 42(3.5+75PrKkr)-¢2(P27.5) (2) 


3.5 Pr Kir? = 6(3.5+75Pr+KKr)- CO (3) 
23 


System bandwidth equation for in = 4: 





_ ————— ee Gl 
[K2-16 (3.5+7SP + Kr)-256]*+ 16(3.5 tKKr2+K~ 16 (P +75) )* 


1/2 


Solving equations (1) through (4): 
K = 252 P=5.93 


ie plies Zod 





24 


4, Sensitivity as a Specification. 

Attention is now directed to the question posed during the 
discussion of example 3.3. If more than one design satisfies 
the specifications, which is the best design? The fact that a 
specific solution was not obtained indicates the need for another 
equation which, for the system being developed, implies the need 
for another specification. Since the design specifications have 
already been incorporated into the solution, it follows that the 
additional specification should be general in nature and 
applicable to all systems. 

One choice for a general control system specification is 


Sensitivity which will be defined as 


i _ dP/p Q P i 
Sa = gel = ails 2 ($a) el 


where Sh is read as "the sensitivity of a function P with respect 
to a parameter Q." (This definition is preferred by the author 
because it imparts more physical meaning than most other 
definitions of sensitivity found in the literature, i.e., it states 
that sensitivity is essentially the ratio of the percentage change 
in the function to the percentage change in the parameter.) The 
reason for this choice is the fact that a sensitivity specification 
will force the final design to be realistic. Specifically, it requires 
the consideration of changes in the parameters of the physical 
system which will result because of tolerances, aging and 
environment. On this basis, the best design is the one which not 
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only satisfies the design specifications, but is also able to 
maintain the desired performance in spite of undesired parameter 
variations. 

Root Sensitivity 

Use of equation (4-1) requires definate of the function P. 
Since sensitivity is being introduced as a means of insuring 
compliance with the design specifications, P should ideally 
be a function of all the design specifications. However, in 
this presentation, only the root specification, i.e., * 4 and 
Wo: will be considered. Hence, the discussion will refer to 
the root sensitivity of the system. (Hereafter, root sensitivity 
and pole sensitivity will be used interchangeably.) 

Restricting attention to root sensitivity does not completely 
preclude the determination or discussion of a best design. In 
fact, there are two arguments which partially justify (and might 
even suggest) this simplification. One argument is that root 
specifications are more stringent than either the bandwidth or 
error coefficient specifications. The latter two are frequently 
| specified in terms of inequalities, i.e., greater or smaller than 
a certain minimum or maximum value. Therefore, some tolerance 
can be built into the design by working on the "safe side” of the 
specifications. Another argument is based on the root- 
coefficient relations of the system characteristic equation. If 


the mobility of two of the systems' closed loop poles are 
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restricted by means of a sensitivity constraint, then the mobility 
of the remaining closed loop poles is also restricted. This 
statement stems from the fact that the specified poles have 
corresponding roots which must be factors of the characteristic 
equation. Now, in order for these roots to be perfect factors of 
the characteristic equation, the coefficients of the characteristic 
equation must have remained reasonably constant for small 
parameter changes. This implies that the remaining roots, and 
hence the remaining poles, also remained reasonably constant. 
Since both bandwidth and the error coefficients are partially 
dependent (and, in fact, totally dependent for low order systems 
without zeroes) on the characteristic equation, it follows that 
their sensitivity functions bear an approximate relation to the 
root sensitivity functions. 
Root Sensitivity Equations 

Although the practical importance of sensitivity has long 
been recognized, it has only been in recent years that effort 
was expended in the development of sensitivity as a useful tool. 
The work of Kokotovic and Siljak, (2), is considered to be a 
major advancement towards this goal. In essence, they used 
the definition of sensitivity, equation (4-1), in conjunction with 
the generalized Mitrovic equations, equations (2-6) and (2-7), 
to obtain algebraic expressions for (what the author considers 


to be) the components of root sensitivity. Their results are not 


a7 


only particularly suited for use in the design of control systems 
based on a specified pole pair, but also sere sensitivity 
computations because they permit the exclusive use of real 
variables in calculating the sensitivity of a function containing 
complex variables. | 

Since the Kokotovic and Siljak sensitivity equations will 
be used in subSequent examples, a brief review of their develop- 
ment will now be presented. Differentiating the generalized 
Mitrovic equations, equations (2-6) and (2-7), with respect to 


a parameter © gives 


tk, tee SKanws, 4, (Ff) + § Gg AZKwr ef aap 


K=O 


K 4-2 
+e OO Pee Cr y= a ( 





a Sik aq wor p x CF) +¢ So a wn b CF) (4-3) 
a Soe wry P« (F) a 


where By ( F )and é,-1' a are the derivatives of Dy, ( F ) and 


jl 


p-1(F ) respectively and are related by the general equation 
oS 
$, (¥) = -2/¢ .., (Ff) US Preoa (F) | 
(4-4) 
Numerical values for D, for various values of e are given in 
Appendix II. Manipulation of equations (4-2) and (4-3) and 
application of equation (4-1) gives 
@) 
yey SS a FB, sy = en 
(4=5) 
LU ry pt ek 
fle SQ t € BzS_ = ~ACe2 (4-6) 
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where 


AA,= a K ax tin bx, (F) ez Z Kank 4.7) (4-7) 
B,= ZaKwk bier (F) Bo =z, On wh $(F) (4-8) 


2 
CE SEK bef) C2 = 222 0Kdelp) 4-9 


Solving equations (4-5) and (4-6) for of and S“n gives 








Q 
wr Le bs (4-10) 
Met — ‘ F i 
‘ww mn ~ Da” (GQ-1] 
= D 
where 
D — a 8, 
A2 Ez (4-12) 
ae - es, 8, 
DE 4 é. ed (4-14) 
ad Se 


Equations (4-10) and (4-11) are the desired sensitivity design 
equations. (Hereafter, 35 and Son will be referred to as the 
“component sensitivity functions.") At this time, their true 
simplicity is masked by their dependence on the solution of a 
number of auxilliary equations. Suffice to state that the solutions 
of the auxilliary equations become exceedingly simple when the 
location of a pole pair is specified ( F , and w,) and the 
coefficients of the characteristic equation are known (ap). This 


will be clearly demonstrated in subsequent examples. 
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Intrepretation of Sensitivity Values 

Equations (4-10) and (4-11) may be used in two ways. 
Substitution of a given set of parameter values into the right 
hand side of the equations gives a definite sensitivity value 
for a specific set of parameter Sie Altermatively, the 
assignment of a sensitivity value to the component sensitivity 
functions give the parameter relations which must be satisfied in 
order to obtain a system of the specified sensitivity. 

In either case, the final design will be dependent upon the 
accepted/assigned value to sensitivity. It follows that some 
guidelines for reasonable sensitivity values must be established. 
Since the literature is lacking in this area, the following guides 
(which are consistent with the original interpretation of equation 
(4-1)) are proposed by the author: 

A, Signs (+ or -) of sensitivity values indicate sense 
only. Therefore, the minimum value of sensitivity is zero, 
implying that a function is totally insensitive to variations of a 
parameLer. 


Bi s 5 = 1 indicates a'linear variation in P due to changes 





in Q.i.e., a 5% change in © results ina 5% change in P. This 
is considered to be the maximum permissible value of sensitivity; 
if possible, systems Should have sensitivity values much less 


than one. 
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oe sg |>2 indicates amplification of the parameter 
variation. This is generally considered undesireable for a 
control system (but may have application in some detection 
schemes). 

D, As indicated in A above, sensitivity values have vector 
properties. Therefore, the total root sensitivity, sé. is 
appropriately expressed as 

(SBT)? = (8G)? + (Syn)? (4-15) 
Sensitivity as a Specification 
The concept of sensitivity has been developed to the extent 


that it can now be used as an additional system constraint 





subject to analysis by the algebraic design method. The 
remainder of this section will be devoted to the design and 
discussion of specific systems which must satisfy sensitivity, 
as well as other, specifications. 
Example 4.1 

Reconsider example 3.3 in view of a specification requiring 
minimum sensitivity of the specified poles with respect to 
variations in system gain, K. 

Considering the system gain, K, as a parameter, the system 


characteristic equation is 


St + s2 (P+ 5+ KKr) + S(SPHKKTE+K) +KZ=O 


(1) 
The sensitivity equations will now be devloped in detail to 
demonstrate the solution of the auxilliary equations. Applying 


equations (4-7) and (4-8) gives 
3) 


A, = 150043 —/00 a2 


Az = 1008 2-5 V2 4 —750V2 43 
. By S00 V2 a3 (2) 


Bz = 1/00 G2 — 2000 


after substituting ee Wie = 5(2)1/2 _— values for 2, (0. Fs) 
and dy 1 (0. 707) from Appendices I and II. Although specific 
values of P and K, are unknown, inspection of equation (1) in 
conjunction with the results of example 3.3 show that the 
coefficients of the characteristic equation must have the 


following values for all solutions 


ee RO a, 2/00 a YTS es 


“t 
wk gt 


which, when substituted into equation (2), gives 


Bi uO fhe *250V2 B, =S0ov2 B3 =-S00 (3) 
Taking the partial of each coefficient in equation (1) with 
respect to K gives 
ote = GY al +kr2 293 = Kz ae © (4) 
Applying equations (4-9) and substituting, as above, for w, 
and 4,(0.707) gives 
C, = soa C2=SoV2 Kr ~SV2(kKrz2 +1) 
which, for Z = 5 as obtained in example 3.3, reduces to 
fe (a2 C. = Sve) (S) 


Applying equation (4-12) and substituting equation (3) gives 





D=-250x/07 (6) 
Applying equation (4-13) and substituting from equations (3) 


and (5) gives 
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Ly, 
Dg B=. 25.00 (7) 


Applying equation (4-14) and substituting from equations (3) 

and (5) gives \ 
De =-/250/2 (1-10Kr) (8) 

Substituting Q=K=50,F = F = 0.707 and equations(6), 


(7) and (8) into equations (4-10) and (4-11) gives 


a = /2 (soKr-7) (9) 


Se7= V2 (10) 

Equations (9) and (10) are the desired design equations. 
Since only one additional constraint equation was needed, they 
would normally be squared and added to obtain an expression 
for the total root sensitivity. In this case, however, SWnis 
observed to be a constant. Therefore, the minimum pole 
sensitivity with respect to K can be found by setting equation 
(10) equal to zero. This gives K, = 0.1 which corresponds to 
the solution of tachometer feedback only. 

Bomk.= 0.1, spn = 1/2 and S, = 0. Therefore, the analysis 
predicts that a ten per cent change in K should produce a five 
per cent change in Wy while F remains constant. These 
predictions can be mathematically verified. The system 
characteristic equation for the specified value of K, is 


Sty s(KKr+S5)r K=O 


For K = 50: Wo = 7.07 F 9 = 9.707 


For K = 55: Ww 7 iA? FF = 0.708 


33 


Therefore Aw» 
Wo 





= 4.95 % zed = 0.1/4 Yo 


Parameter Selection 

The preceeding example raises several questions which 
have yet to be discussed. The first deals with the selection of 
a parameter on which to base the sensitivity analysis. One 
choice might - the parameter which is considered to have the 
widest variations during expected conditions of operation, e.g., 
the gain of an electronic amplifier. Another choice might be 
Similar to the first but restricted to plant parameters. This 
assumes the designer is able to select and/or control the 
sensitivity of the required compensating elements. There is no 
guarantee that either is correct. It may be possible that the 
system is naturally insensitive to the parameter likely to have 
the widest variations in value while it could be highly sensitive 
to a parameter that is expected to have only minute variations 
in value. It is also possible that, by minimizing the sensitivity 
due to one parameter, the effect of another is:amplified., The 
answer, then, requires investigation of all possible parameters, 
It is believed that a thorough investigation of a particular type 
of system could lead to the establishment of general guides for 
the selection of the appropriate parameter(s) in other similar 
systems. Since such an investigation is not consistent with the 


purpose of this paper, the author will be content to arbitrarily 
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select any parameter which will facilitate demonstration of the 
algebraic method. 

It should also be noted that there are no restrictions as to 
the number of parameters which may be similtaneously considered. 
If the sensitivity of pole P with respect to parameters Q), Qo 
and Q2 is desired, the component sensitivity functions for 
each parameter may be individually calculated. The total 
component sensitivity furction due to all three parameters is 


the algebraic sum of the individual components, e.g., 


Wr Wy, Wy 
= s o> 
. since the sense 


of each component must be considered, there is a possibility 
that individual components may cancel. Another fortunate 
circumstance is that the computations for a situation considering 
the root sensitivity with respect to N parameters do not increase 
N-fold. This is because equations (4-7), (4-8) and (4-12), i.e., 
A;, Ag, By, By and D, remain the same for all parameters, 
Minimum Sensitivity 

Another point of interest in example 4.1 is the ease with 
which the condition of minimum sensitivity was obtained. This 
is not a unique situation because it was possible to obtain 
results by inspection in other systems studied also. However, 
it seems quite likely that for many systems the parameter 


relations for minimum sensitivity must be obtained by setting 
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the derivative of the sensitivity function equal to zero. (Note 
the problem if the sensitivity with respect to more than one 
parameter is being considered.) In such cases, care must be 
taken to insure that the function to be differentiated is ina 
form which will indicate minima near zero, This is necessary 
because minimum sensitivity was defined to be zero. The above 
Stipulation can be conveniently accomplished by considering 
either the absolute value or the square of the component 
Sensitivity iynetionmic.., oe or (sgr)* vice sa” 
Use of the total sensitivity function, equation (4-15), is also 
appropriate. 

Discussions thus far have repeatedly used the term "minimum 
sensitivity." This does not necessarily mean that minimum 
sensitivity must be attained by, or even an objective of, every 


design. The specified sensitivity value, whether it be zero or 





two, Should be consistent with the other specifications. For 
example, if the pole locations are only required to be accurate 
within five per cent for satisfactory system response and the 
parameter of interest is also expected to have five per cent 
variations, a sensitivity value of one is sufficient. Of course, 
if analysis shows that a lower value is attainable, the design 
Should be based on the lower value provided it can be achieved 
without additional cost or degradation of some other aspect of 


system performance. 
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The nature of the sensitivity function in the vicinity of the 
design value is also of concern. If the design value is 
associated with a point on the sensitivity function curve where 
the curve is very steep or irregular, the system is likely to be 
more sensitive than an alternative design based on a (reasonably) 
higher value which corresponds to a smooth portion of the 
sensitivity function curve. 

Example 4.2 

A portion of a control system has the forward transfer 
function, G(s) = K/s(s + 2) (s +4). Dominant closed loop poles 
with QP o = 1/2 and wo = 5 are required for the desired response. 
Overshoot is of prime importance; therefore, \e = should be 
maintained within 1% accuracy for .- anticipated gain variations 
of 5%. Velocity and acceleration signals are required elsewhere 
in the system and are therefore available for use as feedback 
Signals. 

The problem as stated gives 3 specifications (pole pair 
and s}.) and 3 parameters. (Note that an additional specification 
concerning Sin would normally require the insertion of an 
additional parameter whereas a specification regarding the total 
pole sensitivity would not. Also note that the specifications 
permit w, to vary, if necessary, in order to maintain - ‘ 


approximately constant.) 
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System CE: s%7 + s2(é+KkKKa)+s(&rKKr)t+K=0 
Mitrovic equations: K= 25 + 25 KKa it) 
K Kr = 22¢+S5 K Ka 


(2) 


Applying equations (4-7) thru (4-9): 


A, = 375 -S50a2 Ai2= S0az-Sa, 
8B, = 250 B2=S0a~-S500 (3) 
Cy = 1-2E GS Cy = S(SKka-kKr) 


Unlike example 4.1, the coefficients of the characteristic 
ondtitiad do not have fixed numerical values. Therefore, the 
solution of the determinants for D and Dy. equations (4-12) and 
(4-14), will be hampered by the presence of all 3 parameters. 


At this point it will be convenient to solve equations (1) and (2) 


for K, and K, - Then equation (3) can be expressed in terms of K 


only. 
. - ke Vo ee ee 
Co esc my 5K (4) 

Therefore 

A), =a 2S ake Go = Koaree 

eo = 250 Be =2K—-Z50 (5) 

C, = 25/K Cz=—10/K 
Using equation (5) to solve for D and Oh: 

é 2 i : (6) 
D =- 4K"+ 500 K- 250 D.= (195 K-13Sx125) 


Using equation (6) to solve equation (4-13: 
sf _ __ 925 (K~ 86.5) (7) 
K Kf£—I125K +125 


Inspection of equation (7) shows that sf = 0 for K = 86,0. 5 oinee 





the only requirement is that sf < 0.2, this gain setting should 
be considered as a possible design value. For K= 86.5, 

i = 0.907 and ane 0284, For these values, the third closed 
loop pole is located at ion = -3.4. Therefore, the complex poles 
are not completely dominant. In addition, figure 4-1 shows the 
sensitivity function is very steep at the point of minimum 
sensitivity. Again using figure 4-1, the system will be designed 
on the basis of K= 600. The other design values are: K, = 0.238 
and K, = 0.0384. The complex poles are assured of being 
dominant because the third pole is now located at 90 = -24. It 
is stated without proof that Syn also approaches zero as K 
becomes very large. It is also less sensitive than SF having 


a value < 0.05 at the design gain setting. 
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9. Conclusions 

The purpose of this paper was to develop a design technique 
based on the premise that design can be accomplished by solving 
a set of algebraic simultaneous equations derived from the 
system design specifications. In addition to demonstrating the 
feasibility of this concept, the examples in the preceeding 
sections also showed that the method developed has certain 
desirable characteristics. Specifically, the algebraic method 
is analytically simple, capable of handling any number of 
parameters, rapid, specific and accurate. 

The method is in no way restricted to the exclusive use of 
the design specifications considered in this paper. The 
specifications considered were simply selected to enhance the 
presentation. In general, a system can be designed to satisfy 
any set of independent specifications providing the specifications 
can be expressed as algebraic equations. It is obvious, of 
course, that the solution of the resulting equations becomes 
more difficult as the number of specifications increase. For this 
reason, it is suggested that analysis be normally limited to 
specifications which are considered to be prime, i.e., those 
which are exactly specified as opposed to others stated in terms 
of inequalities or approximations. The results of this analysis 
can then be used to determine if the secondary specifications are 
also satisfied. 
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The inclusion of sensitivity equations into the set of system 
constraint equations was considered especially cignamleant. Here 
again, the use of sensitivity equations is in no way dependent 
upon the assumptions made in Section 4 concerning the relation 
of root sensitivity to system offh sitivity i.e., use of root 
sensitivity as a performance index. This was logically 
demonstrated in example 4.2 where a root sensitivity component 
was of particular concern while the system sensitivity had no 
restrictions whatsoever. Although only root sensitivity was 
discussed, there is no apparent reason why similar restrictions 
cannot be imposed on other design specifications. 

As previously indicated in Section 4, there is a definite need 
for the study of specification/parameter sensitivity relationships 
before sensitivity can become a useful design tool. At the 
present, the sensitivity effects of all possible parameters must 
be evaluated with respect to each specification and to each other 


before a meaningful design can be achieved. Analysis of even 





a simple system on this basis would require a prohibitive amount 
of computation because of the large possible number of parameter/ 


specification combinations available. 


"Seek Simplicity and Distrust It" 


Alfred North Whitehead 
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